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ABSTRACT 

Most of the previous studies on beryllium abundances in metal-poor stars have 
taken different Galactic populations as a whole when investigating the production 
and evolution of Be. In this Letter, we report on the detection of systematic 
differences in [a/H]-v4(Be) relationships between the low- and high-a stars which 
were identified by previous works. We remind that one should be more careful 
in investigating the Galactic evolution of Be with a sample comprising different 
Galactic populations, because such a mixed sample may lead to inaccurate Be- 
Fe/Be-0 relationships. 

Subject headings: stars: abundances — Galaxy: formation — Galaxy: halo 
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1. Introduction 

During the past two decades, astronomers have been working on beryllium (Be) abun- 
dances in metal-poor stars to investigate the production and evolution of Be in the early 
Galaxy. Most of these studies have taken different Galactic populations (halo and thick-disk 
stars, if included) as a whole when investigating the Be-Fe and Be-0 relationships. This 
may be partly due to that the production of Be is thought to be a global process across 
the Galaxy, and there fore Be abundances should show small scatter a t a given metallicity 



( Suzuki fc Yoshiill200ll ). no matter which population the stars belong to. iBoeseaard fc Novicki 



20061) pres ented, for the first time, some evidences of Be abundances dispersions. ISmiijanic et al 
( 120081 ) and iTan et al.l ( 120091 ) also found some Be-rich stars which are obviously above the 



gener al Be-Fe and Be - tre nds. In their Be abundance survey with the largest sample to 
date, ISmiijanic et al.l (120091 hereafter S09) found a marginal statistical detection of a real 



*Based in part on data obtained from the ESO/ST-ECF Science Archive Facility. 
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scatter in the [Fe/H]-A(Be)Q dia gram. They proposed two possible explanations for this: 
one is that the halo and the thick-disk stars have different [Fe/H]-A(Be) relationships; the 
other is that Be abundances dispersions exist at any given metallicity. 

It is now generally believed that the Galaxy was formed through hierarchical merg- 
ing, which means that different components of the Galaxy may have experienced different 
chemical evolution histories. Since a-elements are mainly produced by core collapse super- 
novae, they are closely associated with the star formation history of the Galaxy. In this 
regard, Galactic components with d ifferent ^-elements abund ances patterns may show dif- 
ferent behaviors in Be abundances. iNissen fc Schuster! (119971 ) found evidence of a bimodal 
distribution of [a /Fe] (a refers to the a verage abundance of Mg, Si, Ca, and Ti) for about a 
dozen halo stars. Gratton et al. (bo03al lbl) found that the accretion component of the Galaxy 
tends to have lower [a/Fe] than the dissipative component. S09 also claimed two distinct 
components of the Galactic halo, but in the log(Be/H)-[a/Fe] diag ram, which seems to make 



the division clearer than in the [Fe/H]-[a/Fe] diagram. Recently, INissen &: Schuster! (12010 



hereafter NS10) performed a precise abundance analysis for a sample of 94 thick-disk and 
halo stars in the solar neighborhood. They confirmed that, in the [Fe/H]-[a/Fe] diagram, 
the halo stars separate into two distinct populations, i.e. the low- and high-a halo popula- 
tions. The very high precision of stellar parameters and a-elements abundances make the 
sample of NS10 ideal for studying the possible different Be abundance patterns for different 
Galactic components in a systematic way. In this Letter, we report on the detection of two 
systematically different Be-a relationships for the stars from the sample of NS10. 



2. The sample 



Among the 94 stars in the sample of NS10, high-resolution and high signal-to-noise 
ratio VLT/UVES spectra covering the Be II 3130 A resonance doublet were available from 
the ESO/ST-ECF Science Archive for 40 stars. We downloaded the reduced spectra for 
th ese stars and the n determined their Be abundances using the same technique as described 
in iTan et al.l (120091). In addit ion, anoth er 3 stars have Be abund ances from the literature (2 



from lRich fc Boesgaardll2009l and 1 from lBoesgaard fc Kingill993l ). For these 3 stars, in order 
to keep our analysis consistent, we first derived their equivalent widths of the Be II 3130 A 
resonance lines based on the Be abundances and stellar parameters from the literature, and 
then re-calculated their Be abundances. For all the sample stars, stellar parameters given by 
NS10 were adopted in the determination of Be abundances. The uncertainty of Be abundance 



1 [A/B] = log[Ar(A)/jV(B)U - log[AT(A)/i\r(B)] , A(X) = log[iV(X)/7V(H)] + 12 
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mainly comes from the uncertainties of stellar parameters and pseudo-continuum location. 
NS10 estimated the uncertainties of effective temperature, surface gravity, and metallicity 
to be ±30 K, ±0.05 dex, and ±0.04 dex, respectively, which leads to a total uncertainty of 
about ±0.03 dex in Be abundance. The typical uncertainty of Be abundance introduced by 
pseudo-continuum location was estimated to be ±0.05 dex. In total, the typical internal 
uncertainty of Be abundance is around ±0.06 dex. Abundances and uncertainties for the 
a-elements were taken from NS10 directly. The abundances results for the sample stars are 
given in Table [TJ We noted that, among the 40 stars whose Be abundances were determined 
directly from the VLT/UVES spectra, 33 stars were included in the sample of S09. If the 
Be abundances given by S09 are scaled to the stellar parameters of NS10, then the average 
difference between our results and that of S09 will be 0.03 ± 0.07 dex (ours minus S09's), 
which is in reasonable agreement within uncertainties. 

In summary, we have consistent Be and a-elements abundances for 43 stars, of which 
14 are thick-disk stars, 13 are low-a halo stars, and 16 are high-a halo stars according 
to NS10 (low-/high-a means low/high [a/Fe] at a given [Fe/H]). As mentioned in NS10, 
the classification of thick-disk and halo is based on kinematic criterion which is somehow 
uncertain (if the velocity distribution of the thick-disk is non-Gaussian with an extended 
tail toward high velocities, then some of the high-a halo stars might belong to the thick- 
disk population). Therefore, in this work we classify the sample stars into two categories 
only, i.e. the low- and high-a stars, based purely on their a-elements abundances. In other 
words, the low-a stars in this work correspond to the low- a halo stars in NS10, and the 
high-a stars in this work include both the thick-disk and the high-a halo stars in NS10. The 
reasonableness of classifying the thick-disk and the high-a halo stars into the same category 
can be further supported by the indistinguishab le abundance patt e rns o f more elements 



for the se two populations as shown in NS10 and iNissen fc Schusterl (120111 ). iGratton et al. 



f|2003al Jbh also mentioned that no clear discontinuity between the rotating part of the halo 



(mostly high-a halo stars) and the thick-disk could be found. 



3. Results and discussion 

Figure [1] shows A(Be) as a function of [Fe/H] and of [a/H] for the sample stars. Here 
we use [a/H] as a proxy for [O/H] as it is generally believed that a-elements have the same 
origin as oxygen. Similar approximation has also been adopted by S09. It can be seen 
clearly in Figure QJa) that Be abundances of the low- a stars are systematically lower than 
that of the high-a stars, and the differences are obviously larger than the uncertainties of 
Be abundances. Even so, one may still expect a uniform A (Be) versus [a/H] relationship 
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for the sample stars because the production of Be is thought to be correlated directly with 
oxygen (and thus a-elements) rather than iron. However, as can be seen in Figure [H(b), Be 
abundances of the low- a stars are still systematically lower than that of the high-a stars. 
As the differences in Be abundances in the [a/H]-v4(Be) diagram are not as obvious as that 
in the [Fe/H]-A(Be) diagram, we made a statistical test to confirm whether the division is 
a real feature or just caused by the uncertainties of the analysis. We made linear fits, in the 
common [a/H] range ( — 1.2 < [a/H] < —0.6), to the high-a stars, to the low-a stars, and to 
the whole sampl^§, and then calculated the root mean square (RMS) deviations from the fits 
respectively. The RMS deviations from the fits for the high- and low-a stars are 0.07 and 
0.06 dex, respectively, which are comparable to the Be abundances uncertainty (0.06 dex); 
while the RMS deviation for the whole sample is 0.18 dex, which is obviously larger than the 
error of Be abundances. Therefore, we conclude that systematic differences in [a/H]-v4(Be) 
relationships do exist between the low- and high-a stars in our sample. One may think 
that Be in the low-a stars have been depleted and thus the abundances are lower than that 
of the high-a stars. This is, however, not the case. As shown in Figure [2j Li abundances 
of the low-a stars are in the upper range o f the sample, b ut their [Be/a] values (solar Be 



abundance A(Be) = 1.32 recommended by iLodderj |2010| was adopted in the calculation 



are systematically lower than that of the high-a stars. As the destruction temperature for 
Li is lower than that of Be, if Be is depleted, then Li should have been depleted more than 
Be. The possibility that the gas from which the low-a stars formed had been exposed to 
the cosmic-rays shorter than that of the high-a stars and thus has lower Be abundances can 
also be excluded, because the low-a stars ar e on average 2-3 Gyr younger than the high-a 



stars (Schuster et al. 2011, in preparation). iPasquini et al.l (120051 ) proposed the possibility 



that some stars may form at very large Galactocentric radii with lower cosmic-ray fluxes 
and thus lower production of Be and heavy elements. However, as shown in Figure OH we 
didn't find any gradient in Be abundances for the sample stars. In the Be abundance survey 
with the largest sample to date by S09, no decreasing trend in Be abundance with increasing 
Galactocentric distance was found, either. 

S09 identified 4 pairs of stars with similar atmospheric parameters and metallicities but 
different Be abundances to testify the existence of scatter in Be abundances. We suggest 
that one of the main reasons for the scatter is the composition of the pairs. Among these 4 
pairs of stars, 3 pairs are included in our sample, and each pair is comprised of one low-a 
star and one high-a star. According to the definition of NS10, at a given [Fe/H], low-a stars 
have lower [a/Fe] (and thus [a/H]) than high-a stars. Moreover, at a given [a/H], low [a/Fe] 



2 All through this Letter, G63-26, HD 132475, and HD 126681 are not included in any linear fits; please 
refer to the last paragraph but one of this Letter for more discussions about these three stars. 



-5 - 



stars have less Be than high [ot/Fe] stars as shown in Figure [T](b). Therefore, it is natural 
that low- a stars have lower Be abundances than high-a stars with the same metallicities. 
Also, it can be easily understood that the separation between the low- and high-a stars in 
the [Fe/H]-v4(Be) diagram is clearer than that in the [a/H]-v4(Be) diagram (compare the two 
panels of Figure [1]). S09 noted the impression of two parallel relations in the Be-Fe diagram. 
However, they did not detect the different [a/H]-A(Be) relationships. This may be partly 
due to the relatively large abundance uncertainties of S09 (note that the stellar parameters 
and a-elements abundances adopted by S09 were taken from different works and had not 
been homogenized). 

As suggested by NS10, the high-a stars may form in regions with rapid chemical evo- 
lution, while the low- a stars may originate from dwarf galaxies with lower star formation 
rates. In this case, the energy spectrum for the cosmic-rays may be different in the regions 
where the low- and high-a stars formed, and thus the production rate of Be could be different 
(as the effective spallation cross-sections are dependent on the energies of the cosmic-rays). 
Therefore, systematic differences in Be abundances between the low- and high-a stars can be 
a natural prediction from the explanation of NS10. We note that the bimodal distribution 



ofJa/Fe] observed by NS10 is consistent with the theoretical results of IZolotov et al.l (12009 



20101 ) . which are based on the assumption that the inner halos of the gal axies contain both 



stars accreted from satellite galaxies and stars formed in situ. However, iNissen &: Schuster 



(120 111 ) shows that abundance patterns of the low-a stars do not match that of the present- 
day dwarf galaxies. They suggest that the low-a stars may have been accreted from more 
massive satellite galaxies at early times. 

As shown in Figure QJb), in the common [a/H] range (—1.2 < [a/H] < —0.6), A(Be) 
of the high-a stars increases with [a/H] faster than that of the low-a stars; however, if the 
"metal-rich" end of the high-a stars are included, then the [a/H]-A(Be) relationship for the 
high-a stars will be flatter than that for the low-a stars. Such a break in the [a/H]-A(Be) 
relationship makes it look like that the "metal-rich" end of the high-a stars are following 
the relation defined by the low-a stars. Similar change in slope also seems to appear in 
the [Fe/H]-v4(Be) relationship at [Fe/H] ~ — 1 as can be seen in Figure Ufa). Inspired 
by the predictio n of sh allower slope at higher metallicity from the mass outflow models, 



Boesgaard et al.l (119991 ) investigated the possible change in slope in the Be-Fe relationship. 
But the metallicity range covered by this work is relatively narrow, and our sample is a bit 
small, hence it is premature to make any conclusions. According to NS10, for halo stars with 
[Fe/H] < — 1.4 ([a/H] < —1.1), the low- and high-a populations will merge together in the 
[Fe/H]-[a/Fe] diagram due to the same origins of iron and a-elements. However, these two 
populations could still be distinguished in the [a/H]-v4(Be) diagram only if the environments 
where these two populations formed had different energy spectrum for the cosmic-rays. In 
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this regard, the [a/H]-A(Be) diagram can be a useful tool to investigate the properties of 
the environments where the stars formed. 

We note that there are three stars standing out th e general Be-a (and Be-Fe) trend 
obviously. HD 132475 was first found to be Be-rich by iBoesgaard fc Novickil (120061) and 



was la ter confirmed by ITan et al.l (120091 ) . HD 126681 was found to be Be-rich by iTan et al. 
(120091 ) . G63-26 is a new Be-rich star found in this work; its Be abundance is about 0.8 dex 



higher than that o f norm al stars with similar [a/H]. According to the results of NS10 and 
Nissen fc Schusterl (120111 ). G63-26 shows normal Na, Mg, Ca, Ti, Cr, Mn, and Cu abun- 
dances, but is overabundant in Si and Ni. The Y and Ba abundances of this star are also in 
the upper range of the sample. The Li abundance of G63-26 is determined to be A (Li) = 2.4, 
which is a bit higher than the Spite plateau. Such a peculiar abundance pattern is very sim- 
ilar to that of HD 106038, w h ich sh ows obvious overabundances of Li, Si, Ni, Y, and Ba as 
described in ISmiljanic et al.l (120081 ) . though the amplitudes of abundance deviations from 
normal stars are much smaller than that of HD 106038. This indicates that the gas from 
which G63-26 formed may have experienced the same but less intensive chemical enrichment 
process as H P 106038 (for example , polluted by nucleosynthetic products of Hypernovae as 
suggested by lSmiljanic et al.ll2008l ). 



In summary, we detect systematic differences in [a/H]-A(Be) relationships between the 
low- and high-a stars in the Galaxy. The [a/H]-A(Be) diagram for different Galactic popu- 
lations can help us understand better the formation and chemical evolution of the Galaxy. 
We remind that, using a sample comprising different Galactic populations to investigate the 
Galactic evolution of Be may lead to inaccurate Be-Fe/Be-0 relationships. For example, in 
the whole metallicity range covered by our sample stars, the equations of linear fits to the 
[a/H]-A(Be) relationships for the high-a stars, for the low- a stars, and for the whole sample 
are A(Be) = (1.07 ± 0.06) [a/H] + (1.35 ±0.04), A(Be) = (1.57 ± 0.17) [a/H] + (1.54 ±0.15), 
and A(Be) = (1.43 ± 0.05) [a/H] + (1.53 ± 0.04), respectively (the corresponding equa- 
tions of the [Fe/H]-A(Be) relationships are A(Be) = (0.98 ± 0.05) [Fe/H] + (1.58 ± 0.05), 
A(Be) = (1.22±0.12)[Fe/H] ± (1.43±0.13), and A(Be) = (1.42 ±0.06) [Fe/H] + (1.88±0.06), 
respectively), which are not the same within uncertainties. However, we should keep in mind 
that the sample of this work is relatively small, and the results need to be confirmed with a 
larger sample. Besides, we use the a-elements abundances as a proxy for oxygen abundances, 
while it is better to investigate the Be-0 relationships directly, though O abundances derived 
from different indicators are usually not consistent. Nevertheless, we hope that the results 
presented in this work could stimulate more research on Galactic evolution of Be and related 
fields, both observational and theoretical. 
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Fig. 1. — Be abundance as a function of (a) [Fe/H] and (b) [a/H] for the sample stars. 
Circles are high-a stars and diamonds are low-a stars. Dash-dotted, dashed, and dotted 
lines are the best linear fits for high-a stars, low-a stars, and the whole sample, respectively. 




Fig. 2. — [Be/a] as a function of A (Li) for the sample stars. The symbols have the same 
meanings as in Figure [TJ the symbols with leftward arrows represent the upper limits of Li 
abundances. 
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Fig. 3. — [Be/a] as a function of (a) perigalactic and (b) apogalactic distances for the sample 
stars. The symbols have the same meanings as in Figure [TJ 
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Table 1. Abundances results, orbital parameters, and classifications for the sample stars 
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-0.59 


0.60 


(!) 


< 1.25 


(2) 


5.82 


8.75 


(2) 


high-a 


HD 121004 


-0.70 


-0.38 


1.02 


(!) 


< 1.34 


(2) 


0.59 


9.01 


(2) 


high-a 


HD 126681 


-1.17 


-0.82 


0.90 


(!) 


1.54 


(2) 


6.14 


8.50 


(2) 


high-a 


HD 132475 


-1.49 


-1.11 


0.60 


(!) 


2.23 


(4) 


3.90 


8.90 


(4) 


(high-a) 


HD 148816 


-0.73 


-0.46 


0.76 


(!) 


1.53 


(2) 


0.70 


9.70 


(2) 


high-a 


HD 160693 


-0.49 


-0.30 


0.80 


(5) 


1.20 


(5) 


2.51 


16.15 


(1) 


high-a 


HD 163810 


-1.20 


-0.99 


0.02 


(!) 


1.67 


(2) 


0.94 


21.21 


(2) 


low-a 
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Table 1 — Continued 



Star 


[Fe/H] 


[a/H] 


A(Be) 


Ref. 


A(U) 


Ref. 




^max 


Ref. 


Class a 


HD 175179 


-0.65 


-0.36 


0.92 


(1) 


< 0.87 


(4) 


1.73 


9.83 


(2) 


high-a 


HD 179626 


-1.04 


-0.73 


0.67 


(1) 


1.81 


(2) 


1.92 


11.42 


(2) 


high-a 


HD 189558 


-1.12 


-0.79 


0.70 


(1) 


2.15 


(2) 


2.34 


9.31 


(2) 


high-a 


HD 193901 


-1.09 


-0.93 


0.01 


(1) 


1.92 


(2) 


0.45 


10.59 


(2) 


low- a 


HD 194598 


-1.09 


-0.91 


0.08 


(1) 


2.02 


(2) 


0.97 


8.88 


(2) 


low- a 


HD 199289 


-1.04 


-0.74 


0.57 


(1) 


2.00 


(2) 


4.94 


8.75 


(2) 


high-a 


HD 205650 


-1.17 


-0.87 


0.41 


(1) 


1.70 


(4) 


3.70 


10.20 


(4) 


high-a 


HD 219617 


-1.45 


-1.17 


-0.37 


(1) 


2.25 


(2) 


1.49 


59.03 


(2) 


(low-a) 


HD 222766 


-0.67 


-0.37 


0.97 


(1) 


< 0.10 


(2) 


1.22 


22.74 


(2) 


high-a 


HD 233511 


-1.55 


-1.21 


-0.17 


(3) 


2.19 


(8) 


0.54 


10.17 


(1) 


(high-a) 
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a Parentheses means the classification is uncertain for stars with [Fe/H] < — 1.4 according to 
NS10. 



